Angular momentum and the geometrical gauge of localized photon states 
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Localized photon states have non-zero angular momentum that varies with the non-unique choice 
of a transverse basis and is changed by gauge transformations of the geometric vector potential a. 
The position operator must depend on the choice of gauge, but a complete gauge transformation of 
a physically distinct state has no observable effects. The potential a has a Dirac string singularity 
that is related to an optical vortex of the electric field. 
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The orbital angular momentum of optical beams has recently received considerable attention 0, . The spin and 
orbital angular momentum of a photon, as traditionally defined, cannot be cleanly separated, although states can 
be constructed that have well defined total angular momentum per photon in any specified direction. Without loss 
of generality we choose this to be the ^-direction, and describe the total, spin and orbital angular momentum along 
this direction with quantum numbers j z , s z and l z respectively. For example, polarized paraxial Bessel beams for 
which j z = are superpositions of two beams: one with l z = 1 and s z = — 1, and one with l z = — 1 and s z = 1 £|. 
In this letter we discuss the angular momentum of photons localized in all three spatial dimensions. We show that 
while geometric gauge transformations can change the angular momenta of the basis states and rotate and rescale 
associated singularities and optical vortices, these features can never be completely eliminated. 

In spite of an extensive literature on nonlocalizability of photons, photon states with arbitrarily fast asymptotic 
power-law Q or exponential [f| falloff of energy density have recently been constructed. What has not been analyzed 
in the past is the deviation of these localized states from spherical symmetry and their consequent angular momentum 
content. It has been argued that a converging or diverging one-photon state can never be localized exactly because of 
mathematical limitations imposed by quantum field theory and, for example, the Paley- Wiener theorem yj. However, 
"^f" ■ a momentum-space basis of exactly localized states such as 

o ■ 

£ r , )A (p)=iVy*e— % A . (1) 

can be constructed with contributions from all p, describing a photon that may be incoming or outgoing relative 
to the spacetime point of localization. This orthogonal basis, while probably not realizable as physical photon states, 
is convenient for calculation of the probability amplitude for photon position and, as we shall show here, for the 
' specification of transverse bases in general. 
Qi The basis states are eigenstates of a photon position operator with commuting components. In spite of a long 
history arguing against the existence of such an operator, we found not one, but a whole family of such position 
operators related by geometric gauge transformations, with the gauge potential defining the rotation of the transverse 
basis about p. Details, together with explanations where arguments against their existence fail, are given in Ref. [||. 
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I We show here that the gauge choice determines the angular momenta of the basis. 

Massless particles possess only two spin states, which can be taken as eigenstates of the helicity operator Sp. 
The resulting coupling of spin and momentum means that the position operator, which generates translations in 
momentum space, generally does not commute with S. For particles of spin 1, the components of S are 3x3 matrices 
that generate rotations of the field vectors, and the position operator is therefore not simply iHV, where V is the 
gradient operator in p-space, but rather a 3 x 3 matrix. One such operator is the Pryce photon position operator 
known since 1948, 

L P = h ^Ip Q Vp~ Q l^pxsj, (2) 

where a = \ for fields and — \ for the vector potential, S is the dimensionless spin-1 operator, and I is the unit 
matrix. In our notation, underscore denotes a matrix and bold denotes a 3-component vector. Thus, the bold face 
signifies that r P has x, y and z components, while its underscore means that each of these components is a 3 x 3 
array that operates on the vector field of a first-quantized photon state, expressed as a 3 x 1 array. This notation is 
carefully maintained to prevent confusion between these two vector roles. 

The Cartesian components of r P do not commute and thus cannot define a basis of localized states. A family of 
position operators that do have commuting components is 

r (x) = ihDp a Vp- a ET 1 (3) 
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where D = e -l ±^e~ l ^- e e~*^2. x is the rotation matrix with Euler angles <j>, 9, \ that rotates the lab z axis into p. The 
role of the matrices D_ and D _1 is to decouple the spin and momentum, allowing the gradient operator to operate on 
the momentum dependence of the field while maintaining the transversality condition. A straightforward calculation 
gives @ 

r (x) = r p - ha M p • S (4) 

with 

a (x) (M) = ^U + V X (M). (5) 
psm9 

The polar and azimuthal angles are denoted 9 and </> in momentum space and and if in position space. 

As the basis vectors for the field and hence for the first-quantized photon wave function, we use complex vectors 
e A of definite helicity A — ±1, with components 

ei*>^)=ei°>,<£)expHAx). (6) 



where 



d M +i\4 u )/j2 (7) 



and we add eo =p to complete the triad. The "hat" denotes a unit vector, and fi = —1,0 and 1, label rows of the 

column vector, e_\, and denote components on the complex vectors (x_ — iy) /V%, £ and (x_ + iy) / \/2, respectively, 
which are eigenvectors of S_ z with eigenvalue fi. Here we express the rotation matrix D. in terms of the same components 
Q and note that e A x ^ (6, 4>) = = D*_^ _ A . The general transverse basis vector e A x ^ is rotated relative to e A °^ by 
the Euler angle \ about p, giving just a phase difference in the helicity basis. 

While the phase of the basis vectors depends on the choice of x, the physical fields are obviously independent of how 
we choose to orient the basis vectors around p. Indeed, we can rotate the basis vectors around p by a different angle 
at different momentum-space positions (9, <j>) , and this cannot change the physical field. In this sense, a reorientation 
transformation \ (9, <j>) — > x' (0j </0 is a true local gauge transformation. It is a basic requirement of the covariance 
of the geometric representation. The invariance of the physical field and hence the photon wave function means that 
the coefficients of the field when expanded in the basis receive compensating phase factors . 

The term may be considered an abelian momentum-space vector potential, analogous to the vector potential 
A of electromagnetic theory 0. The position operator in Eq.Q) depends on the gauge of a^ through Vx (9, <j>) , 
similar to the way the kinetic momentum P of a massive charged particle depends on the gauge of the vector 
potential A. The role of the charge of the massive particle is seen to be taken in momentum space by the helicity of 
the photon, and it is relevant to recall here the well-known result that the helicity defines an invariant subspace of 
the Poincare group. The basis vectors are taken as eigenstates of the position operator, and a gauge transformation 
cannot change their eigenvalues. Thus, a gauge transformation in the basis states e A of the helicity subspace A, say 
e A x ' — > e A x = Te A must change the position operator according to the usual gauge rule 

re A x) r'ei x ' } = Tre A x) 

and this gives the transformation r' = TrT -1 . In our case, T is the phase factor T = e -lA ( x ~ x ) and r= 
ihDp a \7p~ a D~ 1 so that r ' = r — AV (%' — x) j which is exactly the dependence we find for r on the gauge transfor- 
mation. 

The field V x a^ x ^ in momentum space corresponds to that of a magnetic monopole at the origin. The potential 
a^ ) has singular "Dirac" strings of flux lines on the ±z axis that supply the flux emanating from the monopole. 
This is most easily seen by integrating along a path encircling the z axis and equating this to the flux passing 
through the area bounded by the path. The singular strings in a^ x ^ represent an essential nonintegrability or path 
dependence that is responsible for the physical manifestation of the gauge potential [l(J| . As illustrated below, gauge 
transformations induced by changes in x </>) can change the strings, but they do not alter the physical results. As 
shown in Ref. [g, the abelian potential a^ is part of a more general nonabelian gauge potential for SO (3). 

The basis states (01 can be used to express either the ideally localized states £Q or more readily realizable states. 
Adlard, Pike and Sarkar for example, constructed single-photon states with arbitrarily fast asymptotic power-law 
falloff of energy density and photodetection rate and Bialynicki-Birula p| obtained converging or diverging localized 
states with an arbitrarily fast exponential falloff. An advantage of these states is that the falloff rate for the vector 
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potential, the fields, and the Landau-Peierls photon wave function [llj are asymptotically all determined by the same 
exponential factor, and this avoids the problem that the fields themselves associated with exactly localized states are 
not localized 0. As we show now, the gauge choice \ 0) affects the angular momenta of the basis states, whether 
applied to these asymptotically localized states or to the exactly localized states of Rcf . 6] . 
The basis defined by 

X (M) = -mtf (8) 
has total z-angular momentum quantum number j z = to with the single-valued gauge potential 

/„t 7 cos 6 — m 

a (x) =0 r-7-- (9) 

psmt) 

The singularities in a^ ' along the ±z axis (9 — 0, n) are thus changed in strength by the factors 1 =F to. For example, 
for m = 1, the singularity along the positive z axis is missing in a> x ' whereas that along the negative z axis carries 
twice the flux. Other choices of \ ($> 4*) can reorient the singularity along some other direction or replace it by 
a nonintegrable (multivalued) sS- x \ A reorientation of the singularity does not produce any new physics, and as 
discussed above, for simplicity we choose a geometric gauge with the singularity on the ±z-axis. (The most general 
choice of \ can give a singularity that is not straight as discussed in the literature on magnetic monopoles [T^ | , perhaps 
with interesting consequences.) Restricting the Euler angle \ to functions given by Eq.JSJ, the basis vectors can be 
expanded in eigenvectors of the usual spin-1 matrix S_ z and L z — —id/d(j) as 



- 1 _ 2 



(cos 9—1) exp [i (to + 1) (/>] 

sin 9 exp (imcj)) (10) 
(cos 9+1) exp [i (to — 1) q 



The top row (/i = —1) gives the projection of the basis state m onto a state with 5 z -eigenvalue —1 and L z 
eigenvalue to + 1 with probability i (cos 9 — l) 2 ; the second row (/i = 0), has the corresponding eigenvalues and m 

with probability \ (sin# p ) 2 , while the third row (p = 1), has eigenvalues 1 and rn — 1 with probability \ (cos 9 + l) 2 . 
Thus, by inspection, it is confirmed that the total angular-momentum eigenvalue of J_ z of the basis state is m. The 
expectation values of S z and L z for the basis state, obtained from the weighted sum, are then cos 9 and — cos 9 + to, 
respectively, showing that its cosine terms exactly cancel, leaving the eigenvalue to of J_ z 

Restrictions on the uncertainty of the angular momentum of a localized state are imposed by the commutation 
relations between the components of r W and J , which were found in Ref. @ to be 

[JpL k ] = ie jklLk - i\ (d^ x) /d Pk ) . (11) 

Note that the position operator does not transform as a simple vector because, through its coupling to the spin, 
a rotation induces a gauge change. For a photon at the origin for which (r_ k ) = 0, the usual relationship between 
uncertainty and the commutator gives 



dSf ] /d Pk 



(12) 



and 

AJ 2 Ar fc >^(|j^f/^|)- (13) 
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When x is given by Eq.||5J) the z-component of S} x ' reduces to S_ z ~ m ^ = mS-p and within a state space of helicity A, 
dSi'^/dpk = 0. Thus the photon can simultaneously have a definite position and z-component of the total angular 
momentum. However, it does not have definite a; or y-components of J, and there is no definite value for the total 
angular momentum. Nothing can be known definitely about the values of S or L separately. This is consistent with 
the expansion 1)10(1 . 

In coordinate space the electric field describing the localized states discussed here can be written as 



A=±l ' 



x exp [i (p ■ v-pct) /h] (14) 
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where g (9) = sin 9 for the localized states considered in [j| , while g (9) = 1 in Q and and, with the gauge choice 



Jx) fa a\ - Jx) 



^^=^(0,0)^-^. (15) 
To transform to coordinate space we use the expansion in spherical harmonics 

oo I 



cxp (ip • r/fi) = 4^ £ £ (i?, F,"* (9, 0) j, (pr/ft) 

i=0 n=-l 

and integrate over </> to obtain 

oo 

E A*,t) = ^ E i'lT""^) (16) 

x /" d (cos 0) F™^* (0, 0) 5 (0) eg (9, 0) (17) 

where the subscript /i implies the corresponding component in the expansion Hl()[) . The position space field components 
vary as exp [i (to — /S) ip] , indicating a z-component of orbital angular momentum equal to hl z = h (m — /x) . Thus the 
position space z-componcnts of spin, orbital and total angular momentum are exactly the same as those in momentum 
space, and all of the specific results discussed above regarding the angular momentum apply in position space. 
The $ dependence can be obtained by expanding the integrand as 

oo 

v/^ ff (0)4*>,O)= £ c,,Yr"{0,Q) (18) 
l— |m— /j- | 

and using the orthogonality of the spherical harmonics with the same \s z \ value. We consider a few examples. If 
to = 1 and g = 1 then co,i = 4/v / 6! 1 1 <5z i i so that the z-component of the field, ~ sin 79. For the counterclockwise 
rotating component of the field, c\ t i = 4/v / 3^o°^,o +2/v / 6^ / i°<5i,i which gives a independent term and a cos$ term. 
The basis in Ref. 5j implies to = and g — sin 9 and gives Co.; = A/VSY^dtfi — lO/^/SY^Si^- In all cases the field 
component vanishes along any axis for which the corresponding component of L has a nonzero value. 

In place of the common orbital angular momentum operator L = - p x iV in momentum space, we should for 
consistency with our position operator use 

L(x) = x p = DLDT 1 . (19) 

The corresponding spin operator is 

S (x) = J - L (x) = (a&> xplp) Sp (20) 



where J= —iHLp x V + S. The basis vectors are eigenvectors of S-p with eigenvalue A, that is S-pe A = Ae A . They 
are also eigenvectors of the position operator with eigenvalue 0, giving r^eg = and thus L^eg = 0. Thus in a 
basis expansion, LW just differentiates the coefficient of eg, giving no contribution due to the basis. The operator 
S^- 1 alone extracts the total angular momentum of the basis vector eg. Thus use of the position operator, 
separates the angular momentum of the basis from that in its coefficient. The gauge-dependent term & { x) p . S in the 
position operator J1J can best be understood in terms of its relationship to the angular momentum of the basis. In 
position space, orbital angular momentum is associated with a component of the Poynting vector in the tp direction 
such that it spirals along the direction of propagation [l| . In momentum space the functional forms of the position 
and momentum operators are exchanged, and an analogous term appears in the position operator, representing a 
spiraling of the field about p. 

The singular string of discussed above is the axis of a vortex. Expression (|10f) makes explicit the angular 
momenta of the basis vectors along the string, and associated clockwise and counterclockwise rotation about it. The 
polar angle 9 — identifies the positive z-axis and the paraxial limit when describing a beam, while 9 — tt identifies 
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the negative z-axis. If m = 0, the whole z-axis is singular, while if to = 1, there is no singularity associated with the 
positive z-axis (l z — 0), but the negative z-axis has l z — 2, that is it has twice the strength or topological charge. The 
singularity has just been moved from the positive to the negative z-axis. The center of the vortex has zero intensity 
due to the sin 9 dependence discussed above. The orbital angular momentum arises from a bright annular ring about 
the axis, as witnessed in the ji ipr/K) sini? dependence of the field (|16|l . and the radius of this ring goes to zero with 
the parameter describing the spatial extent of the localized photon state. 

In summary, each member of a family of position operators with commuting components defines a corresponding 
basis of transverse unit vectors. The choice of basis contributes a term Aa^ x p to the total angular momentum of 
the basis states and affects the associated optical vortex. However, a complete geometric gauge transformation does 
not change the total field describing a physically distinct photon state. The fact that the position operator given by 
Eq.Q) is not unique is a consequence of unavoidable ambiguity in the selection of a transverse basis, implicit in the 
description of any exactly or partially localized state. 
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